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Motivated by recent angle-resolved photoemission spectroscopy (ARPES) observations of a highly 
metallic two-dimensional electron gas (2DEG) at the (001) vacuum- cleaved surface of SrTiOs and 
the subsequent discussion on the possible role of oxygen vacancies for the appearance of such a 
stated, we analyze by means of density functional theory (DFT) the electronic structure of various 
oxygen- deficient SrTiOs surface slabs. We find a significant surface reconstruction after introducing 
oxygen vacancies and we show that the charges resulting from surface-localized oxygen vacancies 
-independently of the oxygen concentration- redistribute in the surface region and deplete rapidly 
within a few layers from the surface suggesting the formation of a 2DEG. We discuss the underlying 
model emerging from such observations. 



I. INTRODUCTION 

In recent years, SrTiOs (STO), which is a wide band- 
gap semiconductor (3.2 eV), has been commonly used 
as a substrate for the epitaxial growth of oxide het- 
erostructures with other perovskite compounds such as 
the band insulator LaA103 (LAO) or the Mott insula- 
tor LaVOs. In 2004 Ohtomo and Hwang^ detected a 
two-dimensional electron gas (2DEG) state at the inter- 
face of a LAO/STO heterostructure. A polar catastrophe 
was suggested to be the cause for this state. Numerous 
experimental and theoretical studies have been carried 
out on this unusual metallic state since then 2 14 and fas- 
cinating features like superconductivity and coexistence 
with ferromagnetism have been attributed to this 2DEG. 
More recently, Santander-Syro et al 1 and Meevasana et 
al^ found that it is not necessary to have an interface 
with a polar interface with other materials in order to 
obtain a 2DEG at the STO surface. By simply frac- 
turing bulk STO and using angle-resolved photoemission 
spectroscopy, these authors studied a vacuum cleaved 
STO (100) surface and found high carrier densities at 
the bare surface. While the authors' findings opened a 
new scheme in oxide electronics, the precise nature of the 
2DEG has remained elusive. Oxygen vacancies as n-type 
dopants and also as a by-product of the cleaving pro- 
cess were suggested to be related with the 2DEG-. In 
fact, in a recent work by Pavlenko et al. 16 , it was found 
that oxygen vacancies at titanate interfaces are respon- 
sible for a two-dimensional magnetic state and a strong 
orbital reconstruction. 

In view of the previous discussion, we performed 
density functional theory calculations for a number of 
SrTiOs slabs with various concentrations of oxygen va- 
cancies and analyzed the origin of the unusual metallic 
state appearing at the surface of the slabs. The paper 
is organized as follows. In the next section we present 
the computational details of our calculations. Section 3 
is dedicated to the discussion of the electronic structure 
for various slabs and oxygen doping concentrations and 



finally in section 4 we summarize our findings and relate 
them with the experimental observations. 



II. COMPUTATIONAL DETAILS 

We have considered stoichiometric STO slabs on a 
(001) surface. Cleaving bulk STO in the [001] direc- 
tion can lead either to SrO or to Ti02 surface termi- 
nations. Guisinger et al. 17 used cross-sectional scanning 
tunneling microscopy on in situ fractured SrTiOs and 
observed atomically smooth terraces consisting of alter- 
nating surface terminations of Ti02 domains and SrO 
domains. Therefore, both kinds of slabs with SrO and 
Ti02 termination have been included in our study. The 
periodic slabs are separated by a vacuum layer of 20 A. 
In our calculations, we use the experimental lattice con- 
stant 3.905 A from the bulk SrTiOs with a cubic P m3m 
structured The oxygen vacancy is introduced as a defect 
by removing one oxygen atom from the surface. Fig. [T] 
gives an example slab illustrated by a 2 x 2 x 6 supercell 
with (a) SrO and (b) Ti02 termination. Here and in the 
following, we abbreviate the more precise "stoichiomet- 
ric slab with oxygen defect in the SrO surface layer" by 
"SrO terminated slab" . 

In order to clarify the electronic properties of surfaces 
with oxygen defects, we have performed calculations both 
for perfect and for oxygen deficient surfaces. Supercells 
with sizes lxlx 12, 2x2x4, 2x2x6, and 3x3x4 
were designed to represent different oxygen vacancy dop- 
ing concentrations SrTiOs-^ with x = 0.0833, 0.0625, 
0.0417, and 0.0278, respectively. Note that the authors 
of Ref. 1 reported very small oxygen deficiency. We there- 
fore restricted ourselves to concentrations x < 0.1. The 
considered models are also expected to help understand 
the distribution of the extra electrons introduced by the 
oxygen vacancy along the in-plane (xy) and out-of-plane 
(z) directions. 

The calculations were performed in the framework of 
density functional theory (DFT). Structure relaxations 
were done by considering the Born-Oppenheimer pro- 
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FIG. 1. Optimized structures of 2 x 2 x 6 SrTiOs slabs with 
oxygen vacancy in (a) the SrO and (b) the TiCb surface layers. 
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cedure as implemented in the Vienna ab initio simu- 
lation package (VASP) 19 i 2Q . The wave functions were 
represented using the projector augmented wave (PAW) 
metho d 21 ! 22 and we considered the LDA+U functional 23 
with U = 5 and J = 0.64 eV for the Ti 3d states 2 ^ in 
all our relaxations. A plane wave cutoff of 500 eV was 
used. The Brillouin zone was sampled with a /c-mesh of: 
(4 x 4 x 2) for the 2 x 2 x 4 slab, (4 x 4 x 1) for the 2 x 2 x 6 
slab, (lxlxl) for the 3x3x4, and a (8 x 8 x 1) for 
the 1 x 1 x 12 slab. As a small surface oxygen defect con- 
centration is not expected to change the bulk lattice con- 
stant, we relaxed only all internal coordinates of the slab 
but not the lattice parameters. Based on the optimized 
structure, the electronic properties were calculated us- 
ing the all-electron full potential local orbital (FPLO) 25 
method and we employed both the generalized gradient 
approximation (GGA) 26 as well as the GGA+U approx- 
imation as exchange-correlation functionals within DFT. 



III. ELECTRONIC STRUCTURE 

Bulk STO is a band insulator. In Fig. [2] we show the 
GGA calculated band structure of bulk STO with Ti 
3d band weights. The crystal field of the O octahedron 



FIG. 2. Band structure with band weights of Ti for bulk STO. 



surrounding Ti splits the Ti 3d bands into energetically 
lower t2 9 and higher e g states. We observe that the low- 
est unoccupied conduction band consists of three empty 
Ti tig bands which are degenerate at the T point. At low 
temperatures (T < 105 K) STO undergoes a structural 
phase transition to a tetragonally distorted structure 2 ^ 
and this degeneracy of the t^ g bands is lifted. Fig. [2] 
shows the Ti 3d band character analyzed in the cartesian 
coordinate system, with an emphasis on the X — T — X 
and Y — T — Y paths that will become important once 
we break the perfect cubic symmetry by introducing a 
surface. Along the X — V — X path, the band structure 
consists of a weakly dispersive heavy-mass Ti 3d yz band 
and a pair of degenerate strongly dispersive light-mass 
Ti 3d xy and Ti 3d xz bands; along Y — T — Y, electrons 
in the Ti 3d xz band are heavy while those in degenerate 
3d xy and Ti 3d yz are light. 

In the following we will present a comparative analy- 
sis of various STO slabs and discuss the changes in the 
electronic properties observed when considering (i) un- 
doped versus doped STO (ii) unrelaxed versus relaxed 
STO slab structures (hi) GGA versus GGA+U calcula- 



tions, (iv) SrO-terminated versus TiO 2- terminated slabs 
and (v) different oxygen- vacancy concentrations. 

In Fig. [3] (a) we present the GGA calculated density 
of states (DOS) for bulk STO with no oxygen vacancies 
compared to relaxed 2x2x6 SrO and Ti02 terminated 
slabs (Figs. [3] (b) and (c)) with an oxygen- vacancy con- 
centration of x — 0.0417. We observe that the extra 
electrons coming from the oxygen vacancy occupy the 
bottom of the conduction band with mostly Ti 3d char- 
acter, and a metallic state is obtained (compare Fig. [3] 
(a) with Fig. [3] (b) and (c)). One should mention here 
that the (relaxed) STO slabs without oxygen deficien- 
cies show, as in the case of bulk STO, a band insulating 
stated 

Essential for a reliable description of these systems is 
the structural relaxation of the slab geometry. In Fig. 2] 
we show Ti 3d orbital-resolved band structures near the 
Fermi level for unrelaxed and relaxed geometries of the 
1 x 1 x 12 slab with an oxygen- vacancy concentration 
of x = 0.0833 for a Ti02 terminated surface. We ob- 
serve a significant surface reconstruction after structure 
relaxation. While relaxation has little effect on the 3d xz 
band, the splitting between 3d xy and 3d xz is considerable 
after relaxation, and the 3d xz band becomes slightly less 
heavy (m* = 7.9m e (relaxed) against m* = 8.4m e (unre- 
laxed)). This is in contrast to previous calculations that 
were performed without relaxing the structure^. 

In order to analyze the importance of correlation ef- 
fects (at least at the mean-field level as implemented 
in GGA+U) in the oxygen- deficient slabs, we show in 
Fig. [5] the band structure for SrO and Ti02-terminated 
2x2x4 slabs with an oxygen- vacancy concentration of 
x = 0.0625 near the Fermi level region along the T — Y 
path both within GGA and the GGA+U approach with 
U = 5 eV and J = 0.64 eV. Apart from a narrowing of 
the bandwidth in the GGA+U calculations with respect 
to the GGA results (compare Figs. [5] (c) and (d) with 
Figs. [5] (a) and (b) respectively) the main features of the 
metallic state are not influenced. We will concentrate in 
the following on GGA calculations in order to save some 
computational effort but also in order to sidestep an un- 
necessary arbitrariness arising from the choice of U and 
J values. 

Next, we analyze the differences in electronic behav- 
ior of oxygen- deficient SrO-terminated slabs and Ti02- 
terminated slabs. In Fig [6] (a) and (b) we present the 
band structure for a relaxed 2x2x6 slab with an 
oxygen- vacancy concentration of x = 0.0417 for the two 
terminations. Since the slab breaks some symmetries 
present in the bulk, the electronic structure splits into 
a ladder of subbands caused by the Ti atoms becom- 
ing non-equivalent. For both terminations, the Ti bands 
that were unoccupied conduction bands in the bulk (see 
Fig. [2]), are now partially filled with the extra electrons 
gained from the oxygen vacancies, as mentioned above. 
The lowest sub-bands come from in-plane Ti 3d xy states 
largely localized on the surface Ti02 layer (Ti02 ter- 
mination) or on the layer right below the surface layer 
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FIG. 3. Total and partial density of states for bulk STO and 
cleaved slab with oxygen vacancy. (a) bulk STO; (b) 2 x 2 x 6 
with SrO termination; (c) 2 x 2 x 6 with Ti02 termination. 
The insert figures are the zoom in of the bottom of conduction 
band. 



(SrO termination). The slab termination affects to some 
extent the band splitting in the 3d xy bands at T. For 
the same oxygen- vacancy doping, TiO 2 terminated slabs 
lead to the occupation of more Ti 3d xy bands than the 
SrO terminated slabs; this is simply due to the number 
of Ti02 layers significantly affected by the defect which 
is three in case of Ti02 termination versus two in case 
of SrO termination. In fact, a closer look at the slab 
structures reveals that the oxygen- vacancy environment 
in Ti02-terminated slabs is more distorted than in SrO- 
terminated slabs. This trend can also be observed at dif- 
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FIG. 4. Band structure obtained within GGA for a 1 x 1 x 
12 slab (x = 0.0833) with TiC>2 termination, (a) unrelaxed 
structure (b) relaxed structure. 
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FIG. 5. Band structure obtained within GGA and GGA+U 
for a 2 x 2 x 4 slab (x = 0.0625). (a) GGA for SrO termi- 
nation; (b) GGA for Ti0 2 termination; (c) GGA+U for SrO 
termination; (d) GGA+U for Ti02 termination. 
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FIG. 6. Calculated band weights of Ti 3d for the 2 x 2 x 6 
slab with oxygen vacancy, (a) SrO termination; (b) Ti02 
termination. 



ferent oxygen- vacancy concentrations in Figs. [71 (a)-(e) 
where the band structure is shown for relaxed 2x2x4, 
2x2x6 and 3x3x4 slabs (x = 0.0625, x = 0.0417 
and x = 0.0278 respectively) for both SrO and Ti02 ter- 
minated surfaces. Note that the 3x3x4 slabs are com- 
putationally very demanding so that we restricted the 
analysis to the SrO termination. Also, while it appears 
that the Ti 3d xy parabolas are wider for the 3x3x4 
slab (Fig. [7] (e)) compared to the 2 x 2 x 4 and 2x2x6 
slabs (Fig. [7] (a)- (d)), this is only an effect of the different 
Brillouin zones; while the Y = (0, x /2, 0) point is at a dis- 
tance of 0.40 A -1 from the T point in the 2x2 supercells, 
it is at a distance of 0.27 A -1 for the 3x3 supercell. 

Comparing these band structures to ARPES observa- 
tionsii 5 -, we find good agreement. While Meevasana et 
al. see only the light Ti 3d xy bands 15 , Santander-Syro et 
al detect also the shallow heavy bands deriving from Ti 
3d xz . 

We have also investigated the role of the higher lying 
e g states^ after introducing oxygen vacancies and ob- 
serve that in contrast to the LAO/STO system, when an 
oxygen defect is introduced on the bare STO surface, the 
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(a) 2 x2 x4 slab, SrO termination 



(b) 2x2x4 slab, Ti0 2 termination 
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(d) 2x2x6 slab, Ti0 2 termination 
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FIG. 7. Calculated band structures for different (relaxed) 
slabs 2x2x4, and 2x2x6, and 3x3x4 with SrO termination 
(left panels) and TiCb termination (right panels), (a) 2 x 2 x 
4 slab with SrO termination, (b) 2 x 2 x 4 slab with Ti02 
termination, (c) 2x2x6 slab with SrO termination, (d) 
2x2x6 slab with Ti02 termination, (e) 3 x 3 x 4 with SrO 
termination. 



e g states of the Ti atoms neighboring the oxygen vacancy 



are hardly affected. 



A. Two-dimensional electron gas 

Aiura et al 29 studied the behavior of the metallic 
state of lightly electron-doped bulk SrTiOs with angle- 
integrated ultraviolet photoemission spectroscopy (UPS) 
and angle-resolved photoemission spectroscopy. They 
observed a metallic state with a sharp Fermi cut-off in the 
bulk band gap by doping electron carriers. By slight oxy- 
gen exposure, the UPS spectra showed a rapid decrease 
in the spectral intensity of the metallic state followed by 
a sudden energy shift of the valence band. Therefore, 
these spectral features in the bulk system are believed 
to be related to the oxygen deficiency. The difference 
between the metallic state for electron-doped bulk STO 
and electron-doped STO surfaces is the localization of 
the itinerant electrons on the surface in the latter case, 
as we find in our calculations. 

In the following we analyze the confinement of the 
electronic charge in the oxygen- deficient slabs. Figure [8] 
displays the layer-decomposed electron filling of the gap 
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FIG. 8. DOS filling of gap states for layer decomposed PDOS 
of SrO and Ti02 terminated slabs. 



states close to the Fermi level for 1 x 1 x 12 (x = 0.0833), 
2 x 2 x 4 (x = 0.0625) and 2 x 2 x 6 (x = 0.0417) slabs with 
both SrO and Ti02 termination as a funtion of depth. 
The layers are numbered from the bottom of the slab 
and the oxygen vacancy has been introduced on the top 
layer. 

Analyzing Fig. [8] (a), we observe that the number of 
electrons doped into the system due to the oxygen de- 
fect decreases from the surface to the bottom of the slab, 
and there is a dramatic drop after two layers for SrO ter- 
mination and after only one layer for Ti02 termination. 
This indicates that indeed a two-dimensional electron gas 
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(2DEG) forms at the surface. The charge resulting from 
surface-localized oxygen vacancies forms a narrow dis- 
tribution that is peaked at the surface and represents a 
narrow 2DEG in the 3D crystal. The charge distribution 
of Fig. [8] (a) agrees well with the approximate solution of 
the Poisson equation which was used in Ref. [TH to char- 
acterize the observed 2DEG. The number of electrons ef- 
fectively drops to zero after ten layers (if we count TiO 2 
and SrO layers separately) , which means that the system 
recovers its bulk properties within 5 unit cells below the 
surface. With the exception of the top layer in the SrO 
terminated slabs, excess electrons accumulate mostly on 
TiO 2 layers. Nevertheless, the 1 x 1 x 12 slab has an ex- 
cessively large oxygen-defect density at the surface with 
respect to the reported experimental oxygen deficiency 1 . 
Even though the slab is thick enough to see the depth ef- 
fect, the generated top layer has one oxygen vacany per 
one Sr atom (SrO surface termination) and per one TiO 
strip (Ti02 surface termination). This slab was consid- 
ered here for comparison to the calculations performed 
in Ref. [J. 

In Fig. [8] (b), we show the corresponding layer- 
decomposed electron filling of the gap states for a 2 x 2 x 4 
slab. A similar size slab has also been employed as a 
typical model in several previous works We can see 
that the main features are similar to the 1 x 1 x 12 slab 
case, but the bottom layer still shows some electron dop- 
ing. This indicates that this slab is too thin to study 
the carriers in the system if the criterion that a part of 
the slab should show bulk STO carrier densities is con- 
sidered. Fig. [U (c) shows a similar plot for a 2 x 2 x 6 
slab but in this case bulk properties are recovered at the 
bottom of the slab. 

For all considered slabs, the number of doping electrons 
in the gap states are close to 2.0 independent of the size 
of the slab and type of termination. Except for the ex- 
treme case of the 1 x 1 x 12 slab, the second Ti02 layer in 
the SrO surface termination not only contributes maxi- 
mally to the electron distribution but also accumulates a 
similar number of electrons as the 1st Ti02 layer on the 
Ti02 surface termination slabs, while the 2nd SrO layer 
in the Ti02 surface termination slabs has almost no dop- 
ing electrons. This confirms that titanium oxide layers 
play a crucial role for the conductivity in these systems. 

In Fig. [9] we show the Ti 3d charge occupation for the 
SrO-terminated 3x3x4 slab (x = 0.0278). We observe 
that the Ti atoms directly below and diagonally adjacent 
to the O defect position carry a maximal excess charge 
of 0.12 electrons. The other Ti sites in this layer carry a 
rather similar excess charge of 0.11 electrons, making the 
electron doping in this layer almost homogeneous. The 
next Ti02 layer below shows excess charges between 0.04 
and 0.08 electrons, the third Ti0 2 between 0.02 and 0.05 
electrons. This charge distribution in the Ti 3d shells 
indicates that a description of possible correlation effects 
- which is beyond the scope of the present work - would 
require the consideration of site dependent selfenergies. 

A strong instability towards ferromagnetism has been 




FIG. 9. Ti 3d occupations for the SrO terminated 3x3x4 
slab. 



recently discussed 1 — in the context of oxygen- deficient ti- 
tanate interfaces. The density of states (DOS) for our 
2x2x4, 2x2x6, and 3x3x4 slabs for both surface ter- 
minations (Fig. [TO]) shows a sharp peak at ep caused by 
the surface Ti02 layer (Ti02 termination) or the Ti02 
layer below the SrO layer (SrO termination). This in- 
dicates - according to the Stoner criterion - a strong 
instability towards ferromagnetism also in the case of a 
surface 2DEG. Our data shows that by decreasing the 
doping level, i.e. increasing the slab size, the leading 
edge of the valence band shifts to higher energy, and the 
intensity of the sharp peak near the Fermi level decreases 
respectively. 

In Fig. [11] we present the calculated Fermi surface for 
the considered slabs with both SrO and Ti02 termina- 
tions. The figures show that the Ti 3d yz states have 
strongly elliptical Fermi surfaces and high effective mass 
along x compared to low effective mass along y, and 
analogously for Ti 3d xz with the properties along x and 
y exchanged. Meanwhile, Ti 3d xy states have a circu- 
lar Fermi surface and light effective mass along x and 
y. Note that in the calculations, an oxygen defect in 
the SrO surface layer preserves the rotational symmetry 
while a defect in the TiO 2 surface layer breaks it; this is 
why Figs. [11] (b) and (d) are not the same after rotating 
by 90° and exchanging 3d yz by 3d xz . In reality, surface 
oxygen defects would be expected to be randomly dis- 
tributed rather than forming a regular lattice as in the 
calculation. The circular Ti 3d xy derived Fermi surfaces 
compare well to the ARPES measurements of Meevasana 
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FIG. 10. Total DOS for slab 2 x 2 x 4, 2 x 2 x 6, and 3 x 3 x 4, 
SrO and Ti02 terminated normalized to STO formula units. 
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and Santander-Syro et al. 1 . Santander-Syro et al. 
also extract signatures of the elliptical Ti 3d yz /d xz Fermi 
surface pockets^ as found in our calculations. The ellip- 
tical Fermi surface pockets have also been observed by 
Chang et al. — in Ti02 terminated STO single crystals 
where oxygen vacancies were introduced by annealing in 
ultra high vaccuum. 

Finally, we present in Table [I] the effective orbital- 
resolved electronic masses m* for the various slabs calcu- 
lated from 



1 d 2 



h 2 dk 2 



(i) 



where H is the Planck constant and ek is the /c-resolved 
band energy. We find very good agreement with the 
experimental observations for the light effective masses 
(m* = 0.5 ~ 0.7m e ) 1,15 . The heavy effective masses 
are more difficult to determine experimentally and are 
in the range m* ~ 10m e — 20m e 1 ; our values of m* G 
[7.5m e , 9.3m e ] are in reasonable agreement with Ref. [l| 
and in very good agreement with the value m*/m e = 7.0 
measured in Ref. [34|. 

In Fig. [12] we compare the carrier density of our slabs in 
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FIG. 11. Calculated Fermi surface for (a) SrO terminated 
2x2x4 slab; (b) Ti0 2 terminated 2x2x4 slab; (c) SrO 
terminated 2x2x6 slab; (d) Ti0 2 terminated 2x2x6 slab, 
and (e) SrO terminated 3x3x4 slab. 



terms of number of carriers per unit cell area versus num- 
ber of carriers per unit cell volume. Our two-dimensional 
carrier density compares well with the experimental data 
(n 2 D = 2 x 10 14 cm -2 )iii£. We also observe that the 
number of carriers per area doesn't scale with the oxygen 
vacancy concentration as it happens for the bulk electron 
doping risD- This result strongly indicates that the ob- 
served metallic state doesn't correspond to bulk electron 
doping and the 2DEG on the surface is stable indepen- 
dent of oxygen-vacancy concentration and bulk doping 
density. 
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TABLE I. Calculated orbital-resolved effective masses m* for 
different oxygen- deficient slab sizes and surface terminations. 
The error bar for Ti 3d xy is ±0.01, and for Ti 3d xz is ±1. 

lxlx 12 2x2x42x2x63x3x4 



X (ST0 3 -*) 


0.0833 


0.0625 


0.0417 


m* /m e (d xy , SrO) 


0.46 


0.43 


0.43 


m*/m e (d xy , Ti0 2 ) 


0.42 


0.40 


0.47 


m*/m e (d xz , SrO) 


8.6 


8.3 


7.5 


m* /m e (d xz , Ti0 2 ) 


7.9 


9.3 


8.7 



FIG. 12. Calculated carrier density for the SrO and Ti02 
terminated 2x2x6 slabs; (a) carrier density per surface 
area; (b) carrier density per volume. 



IV. SUMMARY 

By employing density functional theory we have stud- 
ied the electronic properties of various oxygen- deficient 
STO surface slabs SrTiOs-^ for x < 0.1 with special 
emphasis on the role of (i) slab structure relaxation, (ii) 
electron correlation at the mean-field level, (iii) surface 



termination and (iv) oxygen-vacancy concentration. We 
detect a significant surface reconstruction after includ- 
ing oxygen vacancies and we observe the formation of a 
metallic state only after introduction of oxygen vacancies. 
The charge carriers - independently of the oxygen con- 
centration - are strongly localized at the surface and de- 
plete rapidly within a few layers from the surface, which 
is an indication of the formation of a two-dimensional 
electron gas. Our calculated Fermi surfaces, effective 
masses and two-dimensional carrier densities show very 
good agreement with experiment. 

The presented calculations show that the presence of 
oxygen vacancies is essential for the formation of a 2DEG 
at the surface of STO and set a basis for further detailed 
investigations on oxygen- deficient surfaces and interfaces 
beyond density functional theory. 
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Appendix A: Additional electronic structure analysis 

We have investigated the contribution to the band 
structure of the Ti atoms in the surface layers containing 
the oxygen vacancy. In Figures [13] and [14j we show this 
analysis for 2 x 2 x 6 slabs with SrO and Ti02 termina- 
tion, respectively. For the SrO terminated slab, the four 
Ti atoms closest to the surface (Ti 1-Ti 4) almost equally 
contribute to the lowest Ti 3d xy band. In the second 
Ti02 layer, contributions from Ti5-Ti8 are again very 
similar, but this time they are responsible for the second 
lowest Ti 3d xy band and to some extent contribute to 
the Ti 3d xz and 3d yz bands. For the Ti02 terminated 
slab (see Fig. [H]) , the fact that the defect breaks the 
rotational symmetry in the surface Ti02 layer is clearly 
visible: Til and Ti2, which are adjacent to the oxy- 
gen defect, show similar contributions to the band struc- 
ture; these contributions are different from Ti3 and Ti4, 
which didn't lose a bond due to the formation of the de- 
fect. Overall, in the Ti02 termination case, the occupied 
bands are less easily assigned to subsequent Ti02 layers 
than in the case of SrO termination. 
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FIG. 13. Calculated band weights of Ti 3d for the 2 x 2 x 6 slab FIG. 14. Calculated band weights of Ti 3d for the 2 x 2 x 6 slab 
with oxygen vacancy and SrO termination. The contribution with oxygen vacancy and TiC>2 termination. The contribution 
of each Ti in the two topmost T1O2 layers is shown separately. of each Ti in the two topmost TiCb layers is shown separately. 
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